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SUMMARY 


Computational methods have been developed to study the trajectories of beta 
particles (positrons) through a magnetic analysis system as a function of the 
spatial distribution of the radionuclides in the beta source, size, and shape of 
the source collimator, and the strength of the analyzer magnetic field. On the 
basis of these methods, the particle flux, their energy spectrum, and source- 
to-target transit times have been calculated for Na^2 positrons as a function 
of the analyzer magnetic field and the size and location of the target. These 
data will be useful in studies requiring parallel beams of positrons of uniform 
energy — such as measurement of the moisture distribution in composite 
materials. 


INTRODUCTION 

Recent studies (refs. 1 to 4) have shown that positron lifetime in compos- 
ite materials is a linear function of the moisture content of the material up 
to the saturation level. This linear dependence of positron lifetime on mois- 
ture content is now being utilized for determining the moisture distribution 
profile in the polymeric targets. Moisture distribution profile is an impor- 
tant factor in determining the environmental degradation of composites. It is 
believed (ref. 5) that nonuniform moisture distribution leads to higher internal 
fatigue of these materials. Ability to measure moisture profile nondestruc- 
tively is of interest since continual absorption and desorption of moisture 
from composites in service usually leaves the moisture profile highly nonuni- 
form inside the test specimen. 

In the previously completed studies (ref. 4) , moisture profiling has been 
accomplished by using a Na22 positron source in conjunction with thin steel beam 
energy degraders. This process is necessarily gross since the positron spectrum 
is continuous in energy. It would be preferable to use positron beams of well- 
defined energy to determine moisture content at the end of their range in the 
test specimen. One practical technique of doing so would require the use of a 
magnetic analyzer for selecting positrons of predetermined energy to probe the 
moisture content at different depths in the test specimen. This paper describes 
analytical methods for maximizing the number of energy-analyzed positrons arriv- 
ing at the target, while minimizing the spread in their energy and arrival time. 
Use of magnetically analyzed positrons would be of practical interest only if 
the spread in their arrival time at the target is less than the time resolution 
of the lifetime measurement system (refs. 1 and 2) and their energy resolution 
is better than their range definition in the target. The computer programs 
needed to calculate various experimental parameters for optimizing the experi- 
mental methods have been developed by Gerald H. Mall and are included as an 
appendix to this paper. 



SYMBOLS 
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W 
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Yc 
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velocity of light 

distance of center of the magnetic field circle from center of source 

diameter of collimator 

kinetic energy of positron 

Fermi function for positrons of energy W 

length of collimator 

particle's path length to target 

slope of the tangent to particle's path at point of exit from magnetic 
field 

particle's rest mass 
particle's relativistic mass 

number of positrons in the energy range from w to W + dW 

transition matrix element for positron decay 

radius of target 

radius of collimator 

radius of magnetic field region 

radius of particle's circular path in the magnetic field 
velocity of particle 

total positron energy in units of moC^ 

maximum value of W, corresponding to end-point beta-spectrum energy 

coordinates of particle; the origin of the x-y coordinate system lies 
at the center of the source spot, with the x-axis passing through 
the center of the magnetic field region 

distance of center of magnetic field from target 

angle of particle emission from source, relative to x-axis 
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T 


transit time of particle 


Tq universal time constant 

(j) angle subtended by particle's circular path in magnetic field at 

circular path center 


PROBLEM STATEMENT 

A well-collimated beam of Na^^ positrons is allowed to enter a variable 
magnetic field. The purpose of this analysis is to follow the trajectories of 
positrons that escape the collimator and enter the magnetic field. The trajec- 
tory of each positron is followed until it arrives at the target located at dif- 
ferent distances fron the center of the magnetic field set to transmit prese- 
lected positron energies. Targets of sizes ranging in diameter from 1 cm to 
2.54 cm, located at a number of distances from the center of the magnetic field, 
are selected for measuring the positron beam. Since the positron emission from 
the Na^^ source nuclei is spherically symmetrical, the calculations are consid- 
ered in a two-dimensional geometry, with the results generalized for a three- 
dimensional geometry. Fran the known energy spectrum of Na^^ positrons 
(ref. 6) , total transit times from the source to the target are calculated for 
each positron arriving at the target. Also, the energy spectrum of the posi- 
trons received at the target location for various settings of the magnetic 
field is determined to assess the definition of the range of the arriving posi- 
trons. Finally, the trajectories of all the positrons that are transmitted 
through preselected magnetic fields are followed for a distance of 20.0 cm from 
the center of the magnetic field. From these data, a distance at which the 
positron energy dispersion and the transit time spread are the smallest and 
the positron number density is the highest is determined for optimal location 
of the target. 


THEORETICAL CALCULATIONS 

For the purpose of developing an analytical basis for the computational 
methods, it is assumed that the positron source is deposited at the bottom of 
a long narrow cylinder. The open end of the cylinder serves to collimate the 
positron beam entering the magnetic field. The source profile is assumed to 
be one of the three practically attainable configurations illustrated in 
figure 1 . The geometrical dimensions of the source collimator considered are 
shown in figure 2. Figure 3 shows a schematic diagram of the experimental 
system. 


Particle's Trajectory Through Source Collimator 

As depicted in the following sketch, the particle is emitted at point 
(0,yT[) at an angle 9 relative to the x-axis with a velocity v; 
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Collimator 


y 



The path of the particle is described by 
X = vt cos 0 
y = vt sin 6 + yi 

Solving equation (la) for the product vt corresponding to x 


(la) 

(lb) 


Lq yields 


vt v= 


x=Lr 


cos 


( 2 ) 


Substituting equation (2) into equation (lb) leads to 


ys = ylx=Lc = ^c tan 0 + yi (3) 

If yg § Rj, or if yg S -r^,, the particle strikes the wall of the collimator 
and is not considered further. If -R^, < yg < Rc» the intersection of the tra- 
jectory with the circular magnetic field is calculated. 


First Intersection of Particle's Trajectory With Magnetic Field 

In the following sketch, the particle's trajectory is represented by the 
line from (0,yi) to (xi,y-|), and the magnetic field by the circle; 
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The equations of the line and the circle are given by 
y = X tan 9 + y^ 

[x - (Rp + Lc)] 2 + y2 = Rp2 (4b) 

Substituting y from equation (4a) into equation (4b) yields 

[~xi - (Rp + Lc)^ 2 + (x^ tan 9 + yi)^ = Rp^ (5) 

Expanding equation (5) and making the substitution Rp + = D leads to 

XI 2(1 + tan2 0) + xi (2y£ tan 9 - 2D) + (d 2 + y^2 _ Rp2) = q (6) 


Equation (6) may be solved for X] ; 
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where the smaller root was selected by subtracting the square root. Expanding 
equation (7) , collecting terms, and reducing the results leads to 


^1 


(D - Yi tan 0) - sec^ 0 - (y^ + D tan 0)2 


sec2 0 


( 8 ) 


Equation (8) gives the x-coordinate of the desired intersection point and the 
y-coordinate may be obtained by substituting equation (8) into equation (4a) . 


Second Intersection of Particle's Trajectory With Magnetic Field 

Upon entering the magnetic field (assumed to be perpendicular to the plane 
of the paper), the particle experiences forces causing it to move in a circular 
path. 
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The forces experienced by the particle satisfy 


qvH = m*v^/Rp 


( 9 ) 


where q is the charge of the particle, m* is the particle's mass at v, and 
H is the magnitude of the magnetic field. The radius of the particle's circu- 
lar path is then given by 


Rp = m*v/qH 


( 10 ) 


and the tangent to the path at the point (x-|,y|) is the slope of the initial 
trajectory, tan 0. The equations of the two circles are given by 


(X - 

D) 2 + y2 = Rp2 

(11a) 

(X - 

Xp) 2 + (y - yp) 2 = Rp2 

(11b) 


Solving equation (11a) for y and substituting into equation (lib) produces 
the following equation for X 2 ; 


(X2 - Xp)2 + ±^p2 _ (X2 - D)2 - 




Rr 


( 12 ) 


Expanding equation (12) leads to 


^2^ - 2x2Xp + Xp2 + Rp2 _ (X2 - D)2 + 2yp^p2 _ (^2 - D) 2 + yp2 = Rp2 


(13) 


By making the substitutions K = 2D - 2xp, L = Xp2 + yp'^ + Rp 
equation (13) can be expressed as 


2 _ p_2 _ n2 


D^, 


KX2 


+ L = ± 2yp^Rp2 - (X2 - D) 2 


(14) 


Squaring both sides of equation (14) and collecting terms leads to the following 
expression for X 2 : 


SX2^ + TX2 + u = 0 


(15) 
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where S = + 4yp^, T = 2 kl - Syp^D, and U = + 4yp^(D^ - Rp2) . Then, 

X 2 is given the following equation where the positive sign of the square root 
has been selected: 


-T + - 


The corresponding value for ^2 determined by substituting equation (16) 

into equation (11a). 


Intersection of Particle's Trajectory With Line x = Lf, + Rp 


The following sketch illustrates the particle's trajectory after leaving 
the magnetic field: 


( 0 , 0 ) 




(X2, Y2)' 


(X3, Y^) 


Magnetic field 



Equation (11b) may be solved for y and differentiated with respect to x to 
find the slope of the tangent to the particle's circular path at (X2,y2)* The 
result is 



m 


X2 - Xp 



(X2 - Xp)2 


0 7) 


The point of intersection of the particle's trajectory with the line 
X = Lq + Rp is then given by 


^3 = ^c + % 

0 8a) 

Y3 = Y2 - ™(^2 - ^c " 

0 8b) 


Determination of Particle's Transit Time 

The particle trajectory along the circular path traverses an angle (|), 
given by 


(j) = cos"^ 


2Rp2 - (X2 - xi)2 - (y 2 - y-i) 


2Rp2 


0 9) 


The total path length is then given by the following expression: 

= Rp(j) + \/(xi - Xi)2 + (y^ - y^)2 + ^(x3 - X2) ^ + (y3 - y2)^ (20) 

For a particle with kinetic energy T, the relativistic mass is given by 
T + mpC^ 

m* = (21) 


where mp is the rest mass and c is the velocity of light. The particle's 
velocity is given by the following equation; 


V 


= \fl - (mo/m*)-^ c 


( 22 ) 


Equations (20) and (22) can be combined to yield the transit time; 


T = Vv 


(23) 
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Using equations developed in these sections the various parameters of 
interest have been calculated for a source, collimator, and magnetic field 
geometry of the type shown in figures 1 and 4, 2, and 3, respectively. The 
results are discussed in the following section. 


DISCUSSION OF TYPICAL RESULTS 

Trajectories of positrons for three different source configurations, trans- 
mitted through the magnetic analysis system set for 200 keV, are shown in fig- 
ures 5 to 7 . It is clear that none of the source configurations considered lead 
to actual positron focusing at distances up to 20 cm from the center of the mag- 
netic field. The flux densities for the flat and convex sources were approxi- 
mately equal and considerably higher than that for the concave source. Also the 
convex source provided slightly better energy and time resolution than the other 
configurations. It was therefore decided to use the convex source configuration 
(fig. 4) for detailed calculations of the transmitted positron beam parameters. 
Seme typical results for the spectra of positrons from a convex source arriving 
at several target locations for two different size targets are shown in figures 8 
to 1 0 and are summarized in tables I to III. These results include the effects 
of Na^2 source deposit profile as well as the positron Coulomb correction factor 
(refs. 7 to 9) . The number of positrons emitted at the source has been taken 
as 1 000 in each case. 

It is clear from the tables that for a 1 -cm-diameter target located a dis- 
tance from the magnetic field center = 12 cm, the dispersion in the arriving 
positron energy is less than 10% of the energy for which the magnetic field is 
set for all cases; that is, Ae/E < 10%. Furthermore, the spread in transit time 
At for positrons arriving at the 1 -cm-diameter target is less than 100 ps for 
Y<, ^ 1 2 cm. This transit time spread is less than the time resolution of the 
lifetime measurement system. For a 2. 54-cm-diameter target, AE/E « 10% and 
At < 1 00 ps for Yj, ^ 20 cm for all magnetic field settings. Thus a target 
to magnetic field center distance of approximately 20 cm appears to be accept- 
able for both the 1-cm and the 2.54-cm targets. This conclusion is also cor- 
roborated by the data shown in figure 11. This figure shows the trajectories 
of all the positrons transmitted to the target through the magnetic analyzer 
when the latter is set for E^+ = 300 keV for 1-cra and 2.54-cm targets. (Sim- 
ilar results are obtained when the magnetic field is set for other positron 
energies . ) 

From the foregoing discussion, it is seen that for the 1-cm target, target 
to magnetic field center distances greater than 1 2 cm would be acceptable, 
though 1 2 cm would yield the highest number density at the target. For the 
2.54-cm target, the minimum acceptable distance is 20 cm. 


CaiCLUSIONS 

A computer program for optimizing a positron magnetic analysis system for 
studying moisture profile in composites has been developed. Three types of 
source configurations were considered. It appears that a convex source deposit 
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provides simultaneously the smallest energy and transit time spread and the larg- 
est positron flux density at the target locations up to 20 cm. For such a 
source, the energy dispersion in the positrons arriving at a 1 -cm-diameter target 
located at or over 1 2 cm from the magnetic field center is less than 1 0% for all 
settings of the analyzing magnetic field. The transit time spread in the posi- 
trons striking the target under these conditions is less than 1 00 ps, which is 
less than the time resolution of the positron lifetime measurement system. An 
examination of the trajectories of the positrons from a convex source transmit- 
ted through the magnetic analyzer for several field settings confirms that a 
separation distance of approximately 1 2 cm also corresponds to the highest posi- 
tron density in the transmitted beam. For a 2.54-cm-diameter target, the opti- 
mum distance would be 20 cm for the smallest energy and transit time spread and 
the largest positron number density at the target. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 20, 1979 
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APPENDIX 


COMPUTER PROGRAMS FOR CALCULATING POSITRON TRANSMISSION 


THROUGH A MAGNETIC ANALYSIS SYSTEM 
Gerald H. Mall 

Conputer Sciences Corporation 
Hampton, Virginia 


PROGRAM FOR CALCULATING ENERGY AND TRANSIT TIMES OF POSITRONS 

Program Description 


The computer program POSTRN is written in FORTRAN IV language for the Con- 
trol Data CYBER 170 series digital computer syston with network operating system 
(NOS) 1.2. The program including interactive graphics routines requires 
60700 octal locations of core storage. A typical case in which 5000 positron 
trajectories are computed requires approximately 20 CPU seconds on the 
CYBER 173. Although the program has been designed to be used interactively, 
it can also be executed in a batch mode. 

The program may be executed in three distinct ways. When the target loca- 
tion, target size, magnetic field setting, and desired number of trajectories 
are specified, the program uses a random number generator and two input proba- 
bility distribution tables to determine the initial conditions of each positron. 
Programs have been developed to generate these tables giving probability versus 
energy and probability versus positron position within the collimator, and these 
programs are described subsequently in this appendix. When the program is exe- 
cuted in this manner, the output consists of several parameters describing those 
positrons which strike the target. Additionally, two disk files (TAPEl and 
TAPE2) are generated which contain sufficient data to reconstruct the trajectory 
of each positron. TAPEl contains data for positrons which miss the target and 
TAPE2 contains data for positrons which hit the target. The second mode of 
operation of the program requires as input one of the disk files TAPEl or TAPE2. 
When this option is selected, program output consists of a plot of the trajec- 
tories of a specified number of positrons which have trajectory data on the disk 
file. The third mode of operation is similar to the second except that the user 
must provide trajectory data for each positron which is to be plotted. These 
data include the magnetic field setting and the energy, position within the 
collimator, and emission angle of each positron. 


Description of FORTRAN Variables 

The following list contains a description of the significant FORTRAN vari- 
ables appearing in the program. The dimension for each array is beside the var- 
iable in parentheses. Each variable is also identified as I, input variable; 

P, program variable; or 0, output variable. Variables which may be input or 
program variables, depending on the mode of operation, are identified by I/P. 
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FORTRAN variable Type 

COLL P 

COLR P 

EFLD I 

EMAX 0 

EMIN 0 

ENER I/P 

MODE I 

NHIT 0 

NMISS O 

NTRY I 

RC I 

RM P 

THETA I/P 

TMAX 0 

TMIN 0 

VEL P 

X(IOOO) P 

Y(IOOO) P 

YC I 

YP I/P 


Description 

Collimator length 
Collimator radius 

Energy for which the magnetic field is set 

Maximum energy of those positrons striking the 
target 

Minimum energy of those positrons striking the 
target 

Positron energy 

Mode of program operation 

Number of positrons which strike the target 
Number of positrons which miss the target 
Number of trajectories to be computed 
Target radius 
Magnetic field radius 
Angle of emitted positron 

Maximum transit time of those positrons striking 
the target 

Minimum transit time of those positrons striking 
target 

Positron velocity 

x-coordinate of points to be plotted 
y-coordinate of points to be plotted 
Distance of target from center of magnetic field 
Position of emitted positron within the collimator 


Flow Charts and Listings of Program 

Program POSTRN .- POSTRN, the main program, performs all input and output 
operations, controls flow through the kinematics and geometry subroutines, and 
calls those subroutines which load the probability tables or generate graphic 
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output depending on the specified mode of program operation. The flow chart 
and listing of POSTRN follow. 
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1 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


PROGRAM PQSTRN(0UTPUT,INPUT»TAPE5“INPUT»TAPE6»0UTPUT> 

1 TAPE1»TAPE2» TAPE3> TAPE4»TAPE7) 

DIMENSION X(1000),Y(1000) 

COMMON X» Y#ENER#XMAS> VEL/H,RM#R P# COLL# COL R» YP# THETA# 

1 XP#X1#Y1»X2»Y2#XCM#YCM#XCP#YCP#X3#Y3#PL#TIME#IJK# 

2 emin#emax#nhit#yc#rc#efld#tmin#tmax#nmiss#mqoe 

DIMENSION EEE(IOO)# YYY( 1 00 ) # TTT ( 100 ) 

COMMON/ JTB/NFR# JR EQ# I BAU D# J CCS# I J 0# T F AC # IJTB ( 4 ) 

NFR * 0 

JREO > 0 

I8AUD =« 120 

JCCS » 0 

TFAC » 0.016 

READ(5#*) MODE 

IF(MODE.NE.O) GO TO 30 

CALL INIT 

READ(5»+) EFLD# YC»RC#NTRY 

UK = 1 

ENER = EFLD 

CALL KINE 

DO 20 I=1#NTRY 

CALL OISTR 

UK » I ♦ 1 

CALL KINE 

CALL GEOM 

IF( X(1 ) .LT.O. ) GO TO 10 
CALL PRT 
GO TO 20 

10 NMISS - NMISS + 1 
20 CONTINUE 

WRITE (6# 1 ) NHIT»NMISS»NTRY#EMIN#EMAX#TMIN#TMAX 
1 FORMAT! IHl# I10»* HITS*#I10#* MI SS ES* # I 1 0# ♦ TRYS*/ 

1 1X#*MINIMUM ENERGY » ♦# E16 .8»5X» ^MAXIMUM ENERGY « *,E16.3/ 

2 1X#^MINIMUM TIME - *# E 16 . 8# 5 X# *MAX I MUM TIME = *#£16. 8) 

STOP 

30 CONTINUE 
CALL PSEUDO 
CALL INIT 

READ(5#*) NTRY,EFLD 
ENER - EFLD 
UK = 1 
CALL KINE 
DO 40 I-1#NTRY 

IF( MODE.EQ.l ) READ(7) E E E ( I ) # YY Y ! I ) # TTT ( I ) 

IF! M0DE.EQ.2 ) READ!5#*) EE E ! I ) # Y Y Y ! I ) # TTT ! I ) 

40 CONTINUE 

DO 50 I-1#NTRY 
UK = I + 1 
ENER = E£E!I) 

YP = YYY!I) 

THETA » TTT!I) 

CALL KINE 
CALL GEOM 

IF!X!l).LT.O. ) GO TO 50 
CALL PLT 
50 CONTINUE 
CALL NFRAME 
CALL CALPLT!0.#0.#999) 

STOP 

END 
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Subroutine INIT .- Subroutine INIT initializes the dimensions of the colli- 
mator and the magnetic field. Several variables used in the output summary are 
also initialized. 



1 


5 


10 


15 


SUBROUTINE INIT 
DIMENSION X( 1000) ,Y{ 1000 ) 

COMMON X>Y,ENER»XMAS»VEL/H»RM^RP#C0LL>C0LR#YP>THETA> 

1 XP#X1» Y1»X2» Y2»XCM» YCM> XC P # YCP, X 3, Y 3# P L » T I ME» I JKf 

2 EMIN# E MAX# NHIT>YC/RC» EFLD# TM IN# T MAX # NM IS S # MODE 
EMIN » 600. 

EMAX - 0. 

TMIN « 100. 

TMAX =• 0. 

NMISS « 0 
NHIT » 0 
RM » 3.48*2.54 
RP = RM 

COLL » 3. *2. 54 
COLR « 0.152*2.54/2. 

RETURN 

END 
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Subroutine KINE .- Subroutine KINE computes the velocity of a relativistic 
positron of known energy. On the first call, the magnetic field strength 
required to focus positrons of the specified energy is determined. On subse- 
quent calls, the radius of each positron's trajectory through the computed 
magnetic field is calculated. 



1 SUBROUTINE KINE 

DIMENSION X{1000),Y(1000) 

COMMON X> Y,ENER, XMAS# VEL»H# RM»RP> C0LL> COL R> YP> THETAt 
1 XP»X1» Y1#X2# Y2#XCM# YCM# XC P # YCP# X3# Y3# P L# T I ME # UK# 

5 2 EMIN# EMAX#NHIT#YC#RC#EFLD# TMIN# TMAX#NMISS# MODE 

T * 1.6E-16*ENER 

XMAS « (T + 9.E+16*9 .llE-31 ) /9.E+16 
TEMP * (9. llE-31/XMAS)+*2 
TEMP - SQRTd. - TEMP) 

10 VEL = TEMP+3.E + 08 

HR « XMAS+VEL/1.6E-19 
HR = l.E+06*HR 
XMAS * 1000. ♦XMAS 
VEL « lOO.+VEL 

15 IF(IJK.EQ.l) H « HR/RP 

IF( IJK.NE.l ) RP « HR/H 

RETURN 

END 
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Subroutine DISTR .- Subroutine DISTR is the executive routine for assigning 
an energy and position within the collimator to each positron subject to user~ 
defined probability distributions. An emission angle is also assigned subject 
to a uniform distribution -tan“^ (2*C0LR/C0LL) ^ 0 ^ tan ^ ( 2 • COLR/COLL) . 



1 SUBROUTINE DISTR 

DIMENSION X(IOOO) »Y( 1000) 

COMMON X, Yf ENER> XMAS>VEL^H> RM»RPf COL L»C 0LR> YP» THETA# 
1 XP#X1#Y1#X2» Y2»XCM# YCM#XCP#YCP# X3# Y3#PL#TIME# UK# 

5 2 EMIN#EMAX#NHIT#YC#RC#EFLD#TMIN#TMAX#NMISS#M0D£ 

CALL EDST 
CALL POSDST 
PI = 0.05061A5 
TMX - -PI 

10 THETA » TMX + 2 . *P I *RANF ( DUM ) 

RETURN 

END 
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Subroutine EDST .- Subroutine EDST assigns a random energy to each positron 
subject to a user-specified probability distribution. On the first entry the 
subroutine reads from TAPES an array of energies and a corresponding array of 
normalized cumulative probabilities. A program which generates these arrays 
is described subsequently in this appendix. A uniform random number between 0 
and 1 is used to interpolate linearly and select an energy. On subsequent 
entries, the read operation is bypassed. 
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1 SUBROUTINE EDST 

DIMENSION X( 1000) »Y( 1000) 

COMMON X# Y/ENERf XMAS# VEL#H/ RM# R P# C DLL# C QL R# YP# THETA# 
1 XP#X1# Y1#X2# Y2#XCH# YCM# XCP # YCP# X 3# Y3# P L# Tl ME # UK# 

5 2 emin#emax#nhit,yc#rc#eflo#tmin#tmax#nmiss#mode 

DIMENSION XX(300)#YY(300 ) 

DATA NIN/0/ 

IF(NIN.NE.O) GO TO 30 
K * 1 

10 NIN * NIN + 1 

REWIND 3 
10 CONTINUE 

READ(3#1) XX(K)#YY(K) 

1 FORMAK 2E16.8) 

15 IF(EOF(3).NE.O) GO TO 20 

K » K + 1 
GO TO 10 
20 K » K - 1 
30 KMl » K - 1 

20 <»0 T - RANF(OUM) 

IF(T.LT.YY( 1) ) GO TO ^0 
DO 50 I»1#KM1 

IF(T.GE.YY{I).AND.T.LE.YY(I+1)) GO TO 60 
50 CONTINUE 

25 I * KMl 

60 ENER - XXCI) + ( XX { I +1 ) -XX ( I ) )* ( T - YY(I))/ 

1 (YY(I+1) - YY( I) ) 

RETURN 

END 



APPENDIX 


Subroutine POSDST .- Subroutine POSDST assigns a random position within the 
collimator to each positron subject to a user-specified probability distribu- 
tion. On the first entry the subroutine reads from TAPE4 an array of positions 
and a corresponding array of normalized cumulative probabilities. A program 
which generates these arrays is described subsequently in this appendix. A uni- 
form random number between 0 and 1 is used to interpolate linearly and select a 
position. On subsequent entries, the read operation is bypassed. 
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SUBROUTINE PQSDST 
DIMENSION X( 1000) »Y( 1000 ) 

COMMON X> Y» ENER#XMAS»VEL,H# RM^ RP> C DLL# C OLR> YP# THETA# 

1 XP#X1» Y1#X2» Y2# XCM» YCM# XC P# YCP# X 3# Y3# P L» TI ME# UK# 

2 EMIN#EMAX»NHIT,YC# RC#EFLD#TMIN#TMAX#NMISS#MODE 
DIMENSION XX(300)#YY(300) 

DATA NIN/0/ 

IF(NIN.NE.O) GO TO 30 
K » 1 

NIN = NIN + 1 
REWIND A 
10 CONTINUE 

READ(4#1) XX(K)#YY(K) 

1 FORMAT( 2E16.8) 

IF(EOF(A),NE.O) GO TO 20 
K « K + 1 
GO TO 10 
20 K - K - 1 
30 KMl = K - 1 
<fO T « RANF(DUM) 

IF(T.LT.YY( 1) ) GO TO AO 
DO 50 I»1#KM1 

IF(T.GE.YY(I).AND.T.LE.YY{I+1)) GO TO 60 
50 CONTINUE 
I » KMl 

60 YP » XX(I) + (XX(I+1)-XX(I) )♦(! - YY(I))/ 

1 (YY{ I+l) - YY( I) ) 

RETURN 

END 
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Subroutine GEOM .- Subroutine GEOM determines the path of each positron. 

For each particle which escapes the collimator, the first point of intersection 
with the magnetic field is calculated. The positron is then assumed to follow 
a circular path until it exits from the field at the second point of intersec- 
tion of the path and magnetic field. From this point the particle travels in 
a straight line until it reaches a distance from the center of the field equal 
to the distance at which the target was placed. The coordinates of this final 
point are calculated as well as the total path length and transit time. 
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SUBROUTINE GEDM 
DIMENSION X( 1000)»Y(1000) 

COMMON X# Y» ENER>XMAS>VEL/H/ RM/ R P/ C OLL> C OLR/ YP# THETA, 

1 XP,X1,Y1,X2,Y2,XCM, YCM, XC P, YCP, X3, Y3, PL, TIME, I JK, 

2 EMIN,EMAX,NHIT,YC,RC,EFLD, TM IN, TM AX ,NM IS S, MODE 
IF(MOOE.NE.O) THETA = TH ETA /5 7. 29577951 
IF(MOOE.NE.O) YC » 20. 

YE » COLL*TAN(THETA) + YP 

X(l) « -1. 

IF( ABS(YE).GT.COLR) RETURN 
X{1) « 0. 

D * RM + COLL 

SECT * 1. + TAN(THETA)**2 

RAO - RM*RM*SECT - ( YP + D*T AN ( THETA ))♦ *2 
XI » (D - YP*TAN(THETA) - SQRT( RAD )) /SECT 
Y1 » X1*TAN(THETA) + YP 
XCM » D 

YCM » 0. 

XCP * XI - RP*SIN(THETA) 

YCP = Y1 + RP*COS(THETA) 

ZK » 2.*(0 - XCP) 

ZL » XCP**2 + YCP**2 + RM**2 - RP**2 - 0**2 

ZS * ZK**2 + <t.*YCP**2 

ZT = 2.*ZK*ZL - 8.*D*YCP**2 

ZU * ZL**2 + A.*YCP**2*( 0**2 - RM**2 ) 

X2 = (-ZT + SQRT(ZT**2 - A . *ZS*ZU ) ) / ( 2 . *ZS ) 

Y2 = SQRT(RM**2 - ( X2 - D)**2) 

DEL * RP**2 - (X2 - XCP)**2 
IFCDEL.LT, 0.00000001) GO TO 10 
SLOPE » (XCP - X2)/SQRT(DEL) 

IFCY2.lt. YCP) SLOPE « -SLOPE 
B » Y2 - SL0PE*X2 
Y3 - YC 

X3 * (Y3 - B)/SLOPE 
GO TO 20 
10 X3 » X2 
Y3 » YC 
20 CONTINUE 

DIST » (X2 - Xl)**2 ♦ (Y2 - Yl)**2 

PHI » ACQS( ( 2.*RP**2 - 0 1 ST ) / (2 . *RP **2 ) ) 

PL * PHI*RP 

PL « PL + SQRT(X1**2 + ( Y1 - YP)**2) 

PL » PL ♦ SQRT((X3 - X2)**2 + (Y3 - Y2)**2) 

TIME » PL/VEL 

RETURN 

END 
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Subroutine PRT .- Subroutine PRT generates two disk files which may be used 
as input to the program when operating in a different mode. If the final coord- 
inates of the positron are within the target area, the energy, initial position 
within the collimator, initial path angle, transit time, and distance from the 
center of the target are output to TAPE2. For those particles which miss the 
target, the information is output to TAPE! . 



1 SUBROUTINE PRT 

DIMENSION X( 1000) #Y( 1000 ) 

COMMON X» it ENER/ XMAS# VEL » H# RM/ R P# C OLL/ C OL R> YP# THETA# 
1 XP#X1#Y1#X2#Y2#XCM# YCM# XC P# YCP# X 3# Y3# PL# TIME# UK# 

5 2 emin#emax»nhit,yc»rc»efld#tmin#tmax#nmiss#mode 

DX =• X3 - COLL - RM 
IF(ABS(DX).LT.RC) GO TO 10 
C A MISS 

THET « THETA+57. 295777951 
10 WRITE(l) ENER#YP»THET#TIME#DX 

NMISS » NMISS + 1 
GO TO 20 
10 CONTINUE 
C A HIT 

15 THET - THETA*57. 295777951 

NHIT « NHIT + 1 
IF( ENER. LT.EMIN) EMIN « ENER 
IF (ENER.GT.EMAX) EMAX» ENER 
IF(TIME.LT.TMIN) TMIN = TIME 
20 IF( TIME.GT.TMAX ) TMAX * TIME 

WRITE{2) ENER# YP#THET#TIME#DX 
20 RETURN 
END 
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Subroutine PLT .- Subroutine PLT produces graphic output. The plot includes 
the collimator, the magnetic field, and the trajectory of each positron. 
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SUBROUTINE PLT 
DIMENSION X( 1000) »Y( 1000) 

COMMON X# Y» ENER# XMAS# VELiH# RM#RP#COLL/C OLR# YP# THETA# 

1 XP#X1#Y1»X2#Y2#XCM» YCM# XC P# YCP# X3# Y3# PL# TIME# UK# 

2 EMIN#EMAX#NHIT#YC# RC#EFLD#TMIN#TMAX#NMISS#MDDE 
DATA JIK/0/ 

JIK » JIK + 1 
IF(JIK.GT.l) GO TO 30 
XMIN ■ 0. 

YMIN * -10. 

XSF « 3. 

YSF » 3. 

X(l) » 0. 

Yd) » -COLR 
X(2) - 0. 

Y(2) - COLR 
X(3) ■ COLL 
Y(3) - COLR 
X(4) » COLL 
Yt4) « -COLR 
X(5) - 0. 

Y(5) » -COLR 
X(6) - XMIN 
Y(6) « YMIN 
X(7) =• XSF 
Y(7) » YSF 

CALL LINPLT{ X, Y# 5# 1# 0#0# 0#0 ) 

XX » COLL 
DX = 0.05 
IK « 2 
X(l) » XX 
Y(l) - 0. 

XF - COLL + 2.*RM 
10 CONTINUE 

XX « XX + DX 

IF(XX.GT.XF) XX » XF 

YY » SQRT(RM**2 - (XX - XCM)**2) 

X(IK) - XX 
Y(IK) =« YY 
IK » IK + 1 
IF(XX.NE.XF) GO TO 10 
XX » XF 
20 CONTINUE 

XX » XX - DX 

IF( XX.LT.COLL) XX » COLL 

YY » -SQRT(RMt*2 - (XX - XCM)*<-2) 

X(IK) * XX 
Y(IK) « YY 
IK » IK + 1 

IF ( XX.NE.COLL ) GO TO 20 
IK = IK - 1 
X(IK+1) » XMIN 
Y(IK+1) » YMIN 
X(IK+2) « XSF 
Y(IK+2) » YSF 

CALL LINPLT(X, Y#IK» 1#0»0#0#0) 
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30 CONTINUE 
X(l) » XP 
Y( 1) « YP 

60 X(2) * XI 

Y(2) > Y1 
IK « 3 
YY - Y1 
^0 CONTINUE 

65 YY • YY + DX 

IF(YY.GT.Y2) YY » Y2 

XX « SQRT(RP**2 - (YY - YCP)**2) + XCP 
X(IK) « XX 

Y(IK) » YY 

70 IK » IK + 1 

IF(YY.NE.Y2) GO TO ^0 
X(IK) « X3 

Y(IK) « Y3 

X(IK+1) = XMIN 
75 Y(IK+1) » YMIN 

X(IK+2) » XSF 

Y(IK+2) « YSF 

CALL LINPLT(X,Y*IK>1»0,0#0# 0) 

RETURN 

80 END 


Other subroutines .- The subroutines PSEUDO, NFRAME, CALPLT, and LINPLT are 
part of the NASA LaRC graphics output library. 


Program Usage 


Input 

All program input is accomplished using FORTRAN list directed reads. The 
data may appear anywhere in the field and when more than one item is specified, 
the data items are separated by commas. The first item expected by the program 
is the variable MODE. Allowed values are 0, 1, and 2. When MODE = 0, the pro- 
gram randomly generates positrons, determines their path through the magnetic 
field, counts the number of particles which strike the target, and, for those 
which do strike the target, determines the minimum and maximum energies and 
transit times. When MODE =1 or MODE = 2 the program produces graphic out- 
put. For MODE = 1, positron trajectories from an earlier MODE = 0 run are 
plotted. For MODE = 2, the user must specify the initial conditions of the 
positrons to be plotted. 

Additional data are a function of the selected mode of operation. When 
MODE =0, the program expects four additional values, separated by commas. The 
first value is the energy, in kilo-electron-volts, for which focusing is 
desired. The program uses this to determine the strength of the magnetic field 
required to focus any positron of this energy on the center of the target. The 
second value is the distance, in centimeters, from the center of the magnetic 
field at which the target is to be placed. The third value is the radius of 
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the target, in centimeters. The final value is the number of positrons to be 
generated. The following is an example of acceptable data for MODE = 0; 


300. ,16. ,0.5, 5000 


For MODE = 1, the user must transfer the data on TAPE! (misses) or TAPE2 
(hits) from a MODE =0 run to TAPE7. The program then expects two additional 
data values. The first value is the number of paths to plot. The second value 
is the energy, in kilo- electron- volts, for which focusing is desired. An 
example follows: 


4,300. 


For MODE = 2, the program expects the same two data items as for MODE = 1 
Additionally, for each path to be plotted the user must provide three items: 
the positron energy, in kilo-electron-volts; the initial y-coordinate of the 
positron within the collimator, in centimeters (-CX)LR ^ YP S COLR) ; and the ini- 
tial path angle, in degrees (-2.9° ^ 9 ^ 2.9°). Sample MODE = 2 input follows 


3,300. 

250. . 0. .-l . 

350.. 0..1. 

300 . . 0 . . 0 . 


Output 

Program output also varies with the mode of operation. For MODE = 0, 
printed output in the following form is generated: 


196 HITS 

MINIMUM ENERGY = 
MINIMUM TIME = 


4804 MISSES 
.28077858E+03 
.12036213E-08 


5000 TRYS 
MAXIMUM ENERGY = 
MAXIMUM TIME = 


.32205860E+03 

.12665825E-08 


The energies are in kilo-electron-volts and the times in seconds. For 
MODE ^ 0 , no printed output is produced, but plots are drawn showing the colli- 
mator, the magnetic field, and the positron paths as illustrated in figure 12. 
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PROGRAMS FOR GENERATING ENERGY AND POSITION DISTRIBUTIONS 

Three programs have been written which generate the probability distribu- 
tion tables required by program POSTRN. Listings of these programs, GENE, GENP, 
and CUMU, are included in this section of the appendix. 

Program GENE tabulates the positron spectrum given by the following expres- 
sion (ref. 10) at intervals of 2.5 keV: 


N(W) dW 



2 

-F(Z,W) (W^ - 


1)V2(Wq - w)2 W dW 


where F(Z,W) for Na^^ has been calculated using the numerical tables given in 
reference 11. Values of F(Z,W) for 25 energy values are given in the GENE 
program listing. The resulting spectrum is illustrated in figure 13. 

Program GENP tabulates the probability of a positron being emitted as a 
function of position within the collimator. It is assumed that the source is 
shaped like a segment of a sphere with a radius equal to the collimator radius 
and a height equal to 25% of the collimator diameter. Results are tabulated 
at 100 equally spaced intervals. The resulting distribution is illustrated in 
figure 14. 

Program CUMU takes as input the output of either program GENE or program 
GENP. Output from CUMU consists of arrays of normalized cumulative probability 
distributions versus energy and position. Output from CUMU is supplied on TAPES 
and TAPE4 as input to program POSTRN. 
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1 PROGRAM GENE(0UTPUT^TAPE6-0UTPUT# INPUT»TAP65«INPUT) 

DIMENSION E(25),F(25)^ X(300)>Y<300) 

COMMDN/JTB/NFR# JREQ>IBAUD> JCCS>rjO»TFAC#IJTB( A) 



NFR - 

0 

5 

JREQ 

■ 0 


IBAUD 

■ 120 


JCCS 

- 0 


TFAC 

- 0.016 


Ed) 

« 10.1 

10 

E(2) 

■ 15.7 


E(3) 

■ 22.5 


E (4) 

■ 30.4 


E(5) 

■ 39.4 


E(6) 

■ 49.3 

15 

E (7) 

- 60.3 


E(8) 

■72.2 


E(9> 

■ 84.9 


E(10) 

- 98.4 


Edl) 

- 112.7 

20 

E(12) 

■ 127.7 


E (13) 

• 143.4 


E(14) 

- 159.6 


E(15) 

- 176.5 


E (16) 

• 193.8 

25 

E(17) 

■ 211.6 


E(18) 

- 246.6 


E (19) 

- 287.2 


E(20) 

■ 327.1 


E(21) 

• 368.1 

30 

E(22) 

- 410.2 


E(23) 

■ 453.1 


E(24) 

■ 496.8 


E (25) 

■ 541.1 


F (1) 

0.2525 

35 

F(2) 

0.3408 


F(3) 

0.4122 


F(4) 

0.4661 


F(5) 

0.5106 


F (6) 

0.5481 

40 

F(7) 

0.5760 


F(6) 

0.6017 


F(9) 

0.6222 


F(10) 

- 0.6414 


F(ll) 

■ 0.6551 

45 

F(12) 

- 0.6684 


F(13) 

■ 0.6797 


F(14) 

• 0.6907 


F(15) 

■ 0.6986 


F(16) 

- 0.7058 

50 

F (17) 

- 0.7130 


F(18) 

- 0.7231 


F(19) 

• 0.7319 


F (20) 

■ 0.7391 


F(21) 

■ 0.7449 

55 

F(22) 

■ 0.7498 


F(23) 

- 0.7539 


F (24) 

■ 0.7574 


F(25) 

■ 0.7602 


WO > 

2.066 
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60 


65 


70 


75 


80 


85 


90 


95 


100 


ENER ■ 12.5 
K ■ 0 
XSUM - 0. 

YSUM - 0. 

10 CONTINUE 

T - ENER/1000. 

W ■ (T + 0.511)/0.511 
XN ■ W**2 - 1. 

XN - W*SQRT(XN)*(WO - W)**2 
DO 20 I-l*2<» 

IF(ENER.6E.E(I).AN0.ENER.LE.E(I-a)) GO TO 30 
20 CONTINUE 
I ■ 24 
30 CONTINUE 

S • F(I) + ((ENER - Ed ))/(E(I+l) - E ( I ) ) ) ♦ ( F ( I +1 ) - F(I)) 
XN ■ S*XN 

WRITE(6d) EN£R,XN 
K - K + 1 
X(K) - ENER 

Y(K) ■ XN*91.538 

XSUM » XSUM + XN 

YSUM ■ YSUM + XN4-ENER 

1 F0RMAT(2E16.8} 

ENER - ENER +2.5 
IF(ENER.LE.540. ) 60 TO 10 
XMIN • 0. 

XSF ■ 60. 

X(K+1) • XMIN 
X(K+2) ■ XSF 
CALL PSEUDO 
Y(K+1) - 0. 

Y(K+2) • 10. 

CALL CALPLT(1.5d.5»-3) 

CALL AXES(0.>0.«0..10.»XMIN» XSF>1 .>0.« 

1 6HENERGY>0.2#-6) 

CALL AXES (0./0.#90.d0.f Y(K+l),Y(K+2)d.»0., 

1 9HFREQUENCY>0.2#9) 

CALL LINPLT(X>Y»K»1»0,0»0,0) 

CALL NFRAME 

CALL CALPLT(0.>0.>999) 

ENAVG • YSUM^XSUM 
URITE(6d) ENAVG 
STOP 
END 
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1 PROGRAM GENP(0UTPUT>TAPE6-QUTPUT, INPUT,TAPE5-INPUT» 

DIMENSION X (300)>Y(300) 

COMMON/JTB/NFR# JREQ>IBAUD> JCCS» IJO>TFAC>IJTB{A) 

NFR • 0 

5 JREO ■ 0 

I6AUD ■ 120 
JCCS ■ 0 
TFAC • 0.016 
K ■ 0 

10 COLR ■ 0.152*2.54/2. 

A ■ COLR 
B ■ A/4. 

BB - (B**2 - A**2)/(2.*B) 

RSQ - A**2 BB**2 

15 DIST ■ -A 

ODIST - A/100. 

10 CONTINUE 

FOIST ■ SQRT(RSQ - DIST**2) + BB 
WRITE(6,1) OIST^FDIST 
20 K - K + 1 

X(K) ■ DIST 
Y(K) • 7.67.6*FDIST 
DIST ■ DIST ♦ DDIST 
IF(DIST.LE.A) GO TO 10 
25 1 F0RMAT(2E16.8) 

XMIN • -0.2 
XSF ■ 0.04 
X(K+1) ■ XMIN 
X(K+2) • XSF 

30 CALL PSEUDO 

Y{K+1) • 0. 

Y(K+2) • 10. 

CALL CALPLT(1.5#1.5*-3) 

CALL AXES(0.>0.>0.»10.f XMIN»XSF»l.f 0.» 

35 1 8HP0SITI0N»0.2>-8) 

CALL AXES (0.,0.,90./10.#Y(K+l)#Y( K+2)»l .#0.> 

1 9HFREQUENCY#0.2,9) 

CALL LINPLT(X> Y*K>1>0*0»0»0) 

CALL NFRAME 

40 CALL CALPLTIO., 0.^999) 

STOP 

END 
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PROGRAM CUMU( OUTPUT# TAP E6-0UTPUT# TAPED 
DIMENSION X(1000)#Y(1000) 

K « 1 

10 CONTINUE 

READ(1#1) X(K)#Y(K) 

1 F0RMAT(2E16.8) 

IF( EQF( 1) .NE.O) GO TO 20 
K ■ K + 1 
GO TO 10 
20 K • K - 1 
DO 30 I-2»K 
Y(I) - Y(I) + YII-l) 

30 CONTINUE 
DO 40 I-1#K 
Yd) - Yd)/Y(K) 

40 CONTINUE 
DO 50 I-1#K 
WRITE(6»1) X(I),Yd) 

50 CONTINUE 
STOP 
END 
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TABLE I.- SUMMARY OF CRITICAL PARAMETERS FOR POSITRONS ARRIVING AT 
TARGET FOR MAGNETIC FIELD SET FOR 200 keV 
Convex source configuration; 1000 positrons onitted at the source] 


Parameter 

Target diameter 




1 cm 

2. 54 cm 


Yj, = 9 cm 


Eg+ (min) , keV 

Eg+ (max) , keV 

(min) , ps 

T (max) , ps 

No. of hits 

175.83 

230.75 

996.10 

1085.64 

68 

1 54.91 
263.40 
982.49 
1118.87 
174 

Yc = 

1 2 cm 


Eg+ (min) , keV 

1 87.34 

165.20 

Eg+ (max) , keV 

218.70 

247.30 

(min) , ps 

1143.53 

1123.25 

T (max) , ps 

121 5.44 

1241 .45 

No. of hits 

44 

140 

K 

o 

II 

1 6 cm 


— 
Eg+ (min) , keV 

1 87.94 

174.86 

Eg+ (max) , keV 

211 .70 

230.75 

(min) , ps 

1350.55 

1320.40 

T (max) , ps 

1411 .32 

1436.76 

No. of hits 

38 

103 

Yc = 

20 cm 


Eg+ (min) , keV 

191 .72 

179.44 

Eg+ (max) , keV 

208.06 

224.71 

(min) , ps 

1549.43 

1518.31 

T (max) , ps 

1587.12 

1610.38 

No. of hits 

28 

77 














TABLE II.- SUMMARY OF CRITICAL PARAMETERS FOR POSITRONS ARRIVING AT 


TARGET FOR MAGNETIC FIELD SET FOR 300 keV 
[convex source configuration; 1000 positrons emitted at the source] 


Parameter 

Target 

diameter 

1 cm 

2.54 cm 


Yc = 

9 cm 


Eq+ (min) , keV 

Eg^. (max) , keV 

(min) , ps 

272.96 

235.58 

347.05 

398.60 

900.02 

894.96 

T (max) , ps 

952.90 

972.04 

No. of hits 

74 

161 


Eg+ (min) , keV 

Eg+ (max) , keV 

(min) , ps 

277.95 

250.92 

328.33 

365.64 

1030.65 

1025.05 

T (max) , ps 

1081 .73 

1101 .24 

No. of hits 

58 

127 


Eg+ (min) , keV 

Eg+ (max) , keV 

(min) , ps 

285.76 

262.36 

317.88 

347.05 

1210.39 

1192.08 

T (max) , ps 

1253.52 

1260.10 

No. of hits 

44 

98 


Eg+ (min) , keV 

Eg^. (max) , keV 

(min) , ps 

288.78 

272.24 

315.12 

332.66 

1381 .06 

1367.96 

T (max) , ps 

1425.30 

1433.93 

No. of hits 

31 

84 



















Target diameter 

Parameter 

1 cm 

2.54 cm 

Yc = 

9 cm 


Eg+ (min) , keV 

Eg+ (max) , keV 

(min) , ps 

352.47 

314.06 

453.85 

511 .48 

851 .16 

846.36 

T (max) , ps 

898.29 

909.75 

No. of hits 

31 

92 


fj, = 1 2 cm 


Ep+ (min), keV 370.30 336.12 

(max), keV 431 .36 472.16 

(min), ps 970.47 970.47 

T (max), ps 1014.94 1023.05 


No. of hits 


Y 



23 64 


1 6 cm 

377.53 

352.47 

414.33 

455.31 

1133.09 

1129.55 

1175.51 

1189.39 

19 

44 


■j, = 20 cm 


383.53 

358.42 

414.33 

434.74 

1293.56 

1288.62 

1336.08 

1342.42 

17 

30 















(a) Convex configuration. 




(c) Concave configuration. 


Figure 


Profiles of Na^2 source deposited at bottom of 
collimator considered in present study. 
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D= 0.386cm 
c 

= 7620 cm 


Figure 2.- Schematic diagram of source collimator. 



Figure 3.- Schematic drawing showing source collimator 
and magnetic field region. 




Figure 4.- Profile of Na^^ source deposited at bottom of 
collimator selected for detailed study. 
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(a) R = 0.50 cm. (b) R = 1.27 cm. 

Figure 10.- Energy spectra of positrons arriving at target located 
at different distances from center of magnetic field set for 
= 400 keV. 
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Figure 14.- Computed probability of a positron being emitted as a 
function of position within collimator. 
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